Scaling of Particle Trajectories on a Lattice I: 

Critical Behavior 



Meng-She Cao and E.G.D. Cohen 
The Rockefeller University 
1230 York Avenue, New York, NY 10021 

February 6, 2008 



Abstract 

The scaling behavior of the closed trajectories of a moving particle gener- 
ated by randomly placed rotators or mirrors on a square or triangular lattice is 
studied numerically For most concentrations of the scatterers the trajectories 
close exponentially fast. For special critical concentrations infinitely extended 
trajectories can occur which exhibit a scaling behavior similar to that of the 
perimeters of percolation clusters. In addition to the two critical exponents 
t = 15/7 and dj = 7/4 found before, the critical exponent o = 3/7, which is 
associated with the scaling function for trajectory size away from criticality, 
also appears. This exponent determines structural scaling properties of closed 
trajectories of finite size when they approach infinity, at criticality. New scal- 
ing behavior was found for the square lattice partially occupied by rotators, 
indicating a different universality class from that of percolation clusters. An 
argument for the scaling behavior found along the critical lines is presented. 
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1 Introduction 



In a number of previous publications^ 1 " 10 ^, the diffusion properties of Lorentz 
lattice gas cellular automata (LLGCA) has been studied. There the behavior of 
a point particle moving through fixed, regularly® or randomly* 1 " 5 ' 7 " 10 ) placed 
scatterers on the lattice sites of a variety of planar lattices has been studied 
numerically. The scatterers consisted either of reflecting right and left mirrors 
or of rotating right and left rotators, which scatter the particle either to its 
right or its left, respectively (Fig. 1-2). The particle is constrained to move 
in unit time steps along the lattice bonds (of unit length) and the lattices 
studied were the square* 1 ' 2 ' 4 ' 7 ^, triangular ( 3 ' 7 \ honeycomb ®, quasi- lattice® 
and random lattice* 3 ' 9 ^. Since these scatterers are fixed on the lattice, the 
particle simply travels periodically in the plane once the trajectory forms a 
closed orbit. In our computer simulations, the particle is stopped right after 
the first period. In almost all cases studied so far, all the trajectories close 
exponentially fast. However, in some cases they close power law slow. While 
in the former case all orbits appear to close after typically 2 10 time steps, 
in the latter case there are extented trajectories, possibly of infinite length, 
which only close after a long (possibly infinite) time. Studying these extented 
trajectories reveals that they exhibit scaling properties, which in the case of 
a lattice fully occupied by scatterers, can be mapped onto a corresponding 
bond or site percolation problem. However, in the case of a lattice not fully 
occupied by scatterers, where empty sites occur and the particle trajectory can 
cross itself, no such mapping seems possible. Nevertheless some of the scaling 
properties of these trajectories are then still the same as those found for the 
fully occupied lattice where no such crossing can occur (Fig. 3). The identity of 
the scaling properties of the closed trajectories with those of a corresponding 
percolation problem is exemplified, for example, by the appearance of the same 
two universal critical exponents that occur in the two dimensional (bond or 
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site) percolation problem: a fractal dimension df = 7/4 associated with the 
length S of the trajectories (perimeter of the percolation cluster) and the mean 
square distance of all points on a large trajectory of length S from the origin, 

R s, 

R 2 S ~ S 2/d f (1) 



as well as the Fisher exponent r = 15/7, characterizing the probability distri- 
bution of closed trajectories of a certain length, 



n s ~ S~ T+i (2) 
A hyperscaling relation holds between these exponents: 

We note that the closed particle trajectories, when mapped onto percola- 
tion clusters, are characterized here by their "surface' ' properties, i.e, their 
length, not, as is usually done, by their "bulk " properties, i.e, the total number 
of lattice sites they contain. The exact values for df were first conjectured to 
be 4/7 by Saproval, Rosso, and Gouyet^ 11 ). Ziff( 12 ), then, developed a scaling 
theory which implied that r should be 15/7 and he verified numerically the val- 
ues of df and r. Later the conjecture was proved by Saleur and Duplantier^ 9 ) 
by mapping the percolation problem onto a Coulomb gas. However, for the 
partially occupied lattice at criticality, there is no theoretical explanation for 
the universality of r and df. 

In this and the following paper, we introduce finer characterizations of the 
closed particle trajectories than we have done in previous publications, where 
only the trajectory size distribution scaling exponent r and the trajectory's 
fractal dimension df were considered. The quantities that describe the finer 
characterizations are: (1) The number of right, left scatterers and empty sites 
on a closed trajectory, i.e. Nr, Nl and Ne, respectively; (2) The winding 



angle W, i.e. the number of right turns minus the number of left turns of the 
particle moving on a closed trajectory; (3) The frequency with which lattice 
sites on the particle trajectory are visited by the moving particle. These are 
more detailed "structural" properties than the "gross" properties incorporated 
in t and df. By studying these properties for finite trajectories of increasing 
length at criticality, we derive, in this paper, a number of trajectory scaling 
properties which include not only the asymptotic behavior for infinitely large 
closed trajectories, but also the approach to the asymptotic behavior. As an 
example of the new properties that we have just introduced, we quote one that 
occurs on the fully occupied square lattice, where each site is occupied either 
by a right or by a left rotator. If Nr and Nl are the numbers of right rotators 
and left rotators contained in a closed trajectory respectively, asymptotically 
Nr and Nl satisfy: 

< ^ - I > c =< ^ - I >^ iV-«- (4) 
N 2 N 2 v ; 

where N = Nr + Nl is the total number of sites of the closed trajectory. 
The averages < ■ • ■ > c and < ■ ■ • > cc are taken over all closed trajectories 
containing N sites in which the particle moves clockwise (Fig. 3) and counter- 
clockwise, respectively. There are two observations to make with regards to 
Eq. (H). First, the number of right and left rotators on a closed orbit become 
asymptotically equal, i.e. asymptotically Nr/N = Nl/N = 1/2. Second 
the critical exponent 0.57 is in good approximation equal to 1 — a, where the 
critical cluster perimeter length exponent a = 3/7. All scaling results observed 
so far for the fully occupied square as well as triangular lattice are of the form 
of Eq.(4), where the constant inside the average (here equal to the critical 
concentration Cr — Cl — 1/2) can vary, while for the partially occupied lattice 
the exponent 0.57 can change to its mean field value 0.50. In the context of 
Eq. (|]), the difference of 0.57 and 0.50 means that asymptotically the right and 
left rotators are almost but not quite randomly distributed over the trajectory, 
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which is not surprising, since some correlations between the placing of right 
and left rotators on the lattice sites to generate a closed trajectory, would be 
expected. Surprising, though, is the universality of the exponent of 0.57, for 
which we have, at present, no explanation. 

Another interesting feature is that some sites on the trajectories are only 
visited once by the moving particle and others can be visited more than once. 
The number of sites visited by the moving particle a different number of times 
exhibits a power law behavior similar to Eq. (||), except for the square lattice 
partially occupied by rotator scatterers, where the scaling behavior differs 
significantly from all other cases. 

The numerical algorithm we used was obtained from an efficient combi- 
nation of the Ziff algorithm^ 4,12,14 ) (use of a virtual lattice) and a technique 
recently developed by Wang and Cohen^ 7 ' 8 -' (dynamical memory allocation). 
The simulation was done on a virtual lattice of size 65536 x 65536. The lattice 
was divided into 1024 x 1024 blocks of 64 x 64 sites rather than 256 x 256 blocks 
of 256 x 256 sites as used by Ziff'- 4 '' and Wang and Cohen^ 7,8 ), since a smaller 
block size is more efficient for dilute scatterer models. 16 bits (2 16 = 65536) 
were used to determine the coordinates (x,y) of a site, the upper 10 bits 
(2 10 = 1024) were used for the location of the block which had been visited 
by the moving particle, while the lower 6 bits (2 6 = 64) were used for the 
site position within the block. Another array of 1024 x 1024 was introduced, 
to record whether a block had been visited by the particle. The application 
of bit shifting, masking etc. to look up the coordinates of blocks and sites 
contained in the blocks made the whole process very fast. Furthermore, we 
used a dynamic memory allocation technique, where an array of pointers is 
generated so that each block has a corresponding pointer and actual memory 
of the states of the sites (i.e. the type of scatterers placed on the site) in the 
block is not assigned to its pointer until the particle enters it (using MALLOC 
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in C). After a trajectory is finished, only memories that have been assigned to 
the pointers are deleted (using FREE in C). 

The advantage of this scheme is that since large trajectories are essentially 
fractals with fractal dimension df = 7/4, only a small fraction of the memory 
for a 65536 x 65536 array is actually used and we do not need to reserve 
memory for those areas that are never visited by the particle. Note that since 
the fluctuations in size of the trajectories are very big, our dynamic memory 
allocation technique saves an enormous amount of memory . For example, 
only a small array is generated for a small trajectory no matter how large 
other trajectories are. After the trajectory is completed, only the blocks that 
have been entered by the trajectories need to be reset to the blank condition 
so that another independent particle can be launched, rather than resetting 
the entire lattice. Therefore the average time required for this operation is 
significantly reduced. Since the memory is allocated dynamically, we do not 
need to know how much memory the largest trajectory requires before we start 
our numerical simulations. 

About 300, 000 independent particles, initially placed randomly on the lat- 
tice, were studied. We only collect closed particle trajectories whose lengths 
are smaller than a certain limit. The trajectory is disregarded if it did not 
close by that number of steps. The value of the limit we used was up to 
2 21 — 2 24 time steps, depending on the concentration of the scatterers. In all 
the simulations, there was no particle that crossed the boundary of the virtual 
lattice. 

We have verified for the fully occupied square lattice that our calculations 
have been made for systems of a sufficiently large number of particles and 
for sufficiently long times that the error bars in our figures are typically of 
the order of the size of the symbols. The same obtains for the fully occupied 
triangular lattice. However, for the partially occupied square and triangular 
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lattices the accuracy of our calculations decreases with the concentration of 
scatterers. 

The organization of this paper is as follows. In section 2, we introduce the 
rotator model on the square lattice and discuss the closed trajectory scaling 
results for the rotator model both for a fully and for a partially occupied 
lattice. In section 3, the same is done for the mirror model. In section 4, we 
discuss the closed trajectory scaling results for the triangular lattice, which 
are the same for both the rotator and the mirror model, since they can be 
mapped in to each other for this lattice ^ 3 ' 7 ^. Section 5 contains a summary 
of our results in three tables as well as a discussion and an outlook. In the 
following paper scaling results, in particular the scaling function, describing 
the scaling behavior in the critical region near criticality will be given. The 
same exponent a that played a dominant role in this paper at criticality, will 
also appear there. 

2 Critical behavior of the rotator model on 
the square lattice 

2.1 The fully occupied lattice 

The rotator model is a special case of a general class of models , originally 
introduced by Gunn and Ortuno* 15 ). Fixed right (Fig. la) and left rotators 
(Fig. lb) are randomly placed on the sites of the lattice. A particle moves 
along the bonds of the lattice and its velocity is rotated either to its right or 
to its left by 7r/2 upon being scattered by a right or a left rotator, respectively 
(Fig. 3a). The total concentration of the scatterers C is the sum of those of 
the right rotators Cr and the left rotators Cl- C = Cr + Cl- 

For the fully occupied lattice, i.e. C — 1, the trajectories of the moving 
particle can be mapped onto the perimeters of bond percolation clusters^ 4 ' 7 - 1 
which was first used by Grassberger* 16 ). A more detailed explanation about 
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this mapping will be presented in the following section by relating the rotator 
model to a mirror model. Since the critical concentration p c for bond perco- 
lation on the square lattice is 1/2, the critical concentration for the rotator 
model is also 1/2, i.e. C Rc = C Lc = 1/2. Note that the right rotators are either 
on the outer side of the trajectories, in which case the trajectory is traversed 
clockwise, or on the inner side of the trajectories, in which case the trajectory 
is traversed counter-clockwise (Fig. 3a). It is convenient, in our numerical 
simulations, to introduce for the analysis of the trajectories the winding angle 
W, which, for example allows a determination of whether a closed trajectory 
is traversed clockwise or counter-clockwise. For a given trajectory, W is com- 
puted, at each step, by increasing it by n/2 if the particle is turned to the 
right, and decreasing it by tt/2 if the particle is turned to the left. Then, 
when the trajectory closes, the winding angle will either be 2n or — 2n which 
corresponds to clockwise or counter-clockwise rotation , respectively, because 
trajectories cannot cross themselves. 

Although we expect that, for large trajectories, the number of right rotators 
Nr and the number of left rotators Nl contained in the trajectories are on 
average the same, there are symmetric fluctuations of Nr/N with respect to 
the mean value Cr c = 1/2. The distribution of N R /N can be fitted to a 
double Gaussian (Fig. 4a). The symmetry is due to the fact that the same 
trajectory can be generated by replacing all right (left) rotators with left (right) 
rotators, respectively, and reversing the particle velocity. However, if we just 
look at the clockwise traversed trajectories, we find that the fluctuations of 
N R /N are no longer symmetric with respect to 1/2 (Fig. 4b). This is not 
difficult to understand, because the right rotators are always on the outer side 
of the trajectories while the left rotators are always on the inner side of the 
trajectories. In other words, the symmetry between right rotators and left 
rotators is broken by our selection. From our numerical results, we found 
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that for most clockwise closed trajectories Nr > Nl, while in a few cases, 
Nr < Nl (Fig. 4b). As we increase the size of the trajectories, the average of 
the distribution of Nr and Nl shifts toward 1/2 with a power law (Fig. 5), 

< ^ - 1/2 > c =< f - 1/2 > cc ~ N~™ (5) 

It has been have noticed before^ 4 ' 9 ) that memory effects play a dominant 
role in the motion of the particle on the lattice with fixed scatterers and in 
the generation of its trajectory. However, no direct measurement of a memory 
effect has been given. In this paper, we do so for the rotator model, by 
examining how many sites on a trajectory are visited once Ni (no memory 
effect) and how many sites on the same trajectory are visited twice N 2 (memory 
effect). A site on the trajectory can not be visited more than twice on the fully 
occupied lattice, since for a particle to return to the same site, it has to move 
an even number of steps. As at each step, the winding angle changes by ±7r/2, 
the change of the winding angle can only be either 2n or n when the particle 
returns to the same site, (Fig. 3), so that only Ni and N 2 are allowed. Thus 
we have the following sum rule, 

Ni + N 2 = N (6) 

If S is the length of the trajectory, we also have, 

Ni + 2N 2 = S (7) 

Our numerical results show that (Fig. 6), 

Ni 1 1 N 2 nw 

< - o >=< o - "77 >~ N 8 
N 2 2 N w 

Note that Ni/N approaches 1/2 from above while N 2 /N approaches 1/2 from 
below, i.e. for infinitely large trajectories Ni/N and N 2 /N are equal. This 
implies that memory effects are indeed very important even for those trajec- 
tories whose length are infinitely large. The asymptotic value 1/2 for Ni/N 
and N 2 /N was noticed independently by Ziff( 17 ). 
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Besides the memory effect exhibited by Ni and N 2 , there are other interest- 
ing structural properties associated with Ni and N 2 . For each site belonging 
to Ni, only two out of four adjacent bonds belong to the trajectory, while for 
each site belonging to N 2 , all four adjacent bonds belong to the trajectory (Fig. 
3a), i.e. the site is surrounded by the trajectory (Fig. 3a). In other words, 
the Ni form the "surface" of a trajectory, while the N 2 form the "bulk" of the 
trajectory. We want to point out that the "surface" and the "bulk" introduced 
here are different from the outer and inner sides we introduced earlier, since 
some sites on the outer side of a trajectory belong to N 2 and some sites on the 
inner side of that trajectory belong to Ni. 

Although Eq. (§) and Eq. @ exhibit the same power law behavior, we 
have not obtained an explanation for this. From Eq. ([]), Eq. ([?[) and Eq. (|8|), 
the following relation can be derived, 

< § - I >~ iV-°' 57 0) 

indicating that S and N are asymptotically proportional to each other. There- 
fore, replacing S~ ' 57 by iV _a57 in all scaling equations is justified, as long as 
the asymptotic behavior is considered. 

2.2 The partially occupied lattice 

For a partially occupied lattice, the trajectories can cross themselves (Fig. 
2a). The outer side and inner side of the closed orbits cannot be well defined 
anymore, so that a direct mapping of the trajectories onto the perimeters of 
bond percolation clusters is not possible^ 7 ' 18 -'. Nevertheless earlier numerical 
studies^ 7 -* have still shown the existence of two critical lines symmetric with 
respect to the line Cr = Cl (Fig. 7) and the critical exponents, r and df 
have the same values as at C = 1. However, the existence of these critical 
lines below C = 0.56 could not be established, due to the prohibitively long 
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numerical calculations needed. We will give more details about some special 
features of the critical lines in the following paper. 

One important consequence of the inequality of Cr and C L at criticality for 
C < 1 is that the winding angles are asymptotically proportional to the number 
of sites or to the number of time steps contained in the closed trajectories. 
As an example, we consider to the critical point at concentration C = 0.90, 
namely, Cr = 0.477, Cl = 0.423. To our surprise, our numerical results 
exhibit that the deviation of the winding angle from its mean value (Cr — Cl) 
decays not with a power law but with a stretched exponential law (Fig. 8), 

<|^-(^-C L )|>^2(- 2 -^ - 18 ) (10) 

For Nr/N, N l /N and N E /N, our numerical calculations yield the following 
scaling relation (Fig. 9), 

^ \ Nr r \ ^ \ Nl n I ^ ^ \ Ne n I ^ j\r-o.5o /-.-n 
< |— - C R \ >~< |— - C L \ >~< |— - C E \ >~ TV (11) 

Note that the value of the exponent —0.50 differs from that for the fully 
occupied lattice —0.57, indicating that self-crossing occurs randomly as the 
particle generates its trajectory. 

The memory effects on the partially occupied lattice are more complicated 
than on the fully occupied lattice. For example, the number of times that a 
site can be visited by the moving particle can range from one to four. We 
denote the number of these different types of sites by N 1: N 2 , N 3 and N 4 , 
respectively. An interesting feature for C < 1 is that "reflectors" can be 
formed after a number of time steps on a trajectory (a "reflector" is formed 
when the moving particle travels through the same bond as it has traveled on 
before, but in the opposite direction , Fig. 10. We remark here that these 
"reflectors" , are responsible for the absence of diffusion for the flipping rotator 
model as proved by Bunimovich and Troubetzkoy^ 19 ^ and make N 3 and A" 4 non 
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vanishing. Since each site belongs to only one of the four types, Ni, N 2 , N 3 
or N4, we have the following two sum rules, 

N 1 + N 2 +N 3 + N± = N (12) 

and 

Ni + 2N 2 + 3N 3 + 4N 4 = S (13) 

These two relations are quite useful for numerical simulations, since for known 
Ni N 2 , N and S, one can find N 3 and A/4, saving therefore a large amount 
memory and computational time. If the asymptotic values for N±/N, N 2 /N, 
N 3 /N and N^/N are represented by K±, K 2 , K 3 and K4, respectively, our 
numerical simulations show the following power law behavior (Fig. 11), 

< — — Ki >~< — - K 2 >~< K 3 - — >~< K A - *± >~ iV- 39 (14) 
N N N N K J 

The sum of Ki, K 2 , K 3 and K<± is equal to one as required by Eq. (|12j). We 
found that Ki, K 2 , K 3 and K4 are functions of the concentration C and all 
converge to values near 1/4, (Fig. 12). The exponent —0.39 in the above equa- 
tion is significantly different from the corresponding one for the fully occupied 
lattice —0.57, although r and df are still the same as found by Cohen and 
Wang (7) . The same scaling behavior is found at C = 0.80, suggesting that the 
critical behavior is universal along the critical line. 

In the same way as was done for the fully occupied lattice, we can derive 
the scaling behavior for S/N from Eq. (|T^), fll^D and (|bj|), 

S 

<—- 2.3670 >~ iV -0 ' 39 (15) 

We note that the sites that belong to N 2 can be separated into two classes: 
those that only have three adjacent bonds traversed by the particle and those 
that have four adjacent bonds traversed by the particle. The former together 
with Ni form the "surface" of the trajectories while the latter together with 
A3 and N4 form the "bulk" of the trajectories. 
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3 Critical behavior of the mirror model on the 
square lattice 

3.1 The fully occupied lattice 

The mirror model was proposed by Ruijgrok and Cohen sometime ago ^\ 
where the scatterers consist of right mirrors (tilted to the right by 7r/4, Fig. 
lc) and left mirrors (tilted to the left by 7r/4, Fig. Id), which reflect the 
particle upon collision, like a photon is reflected from a mirror. Cr and Cl 
are now the concentrations of right mirrors and left mirrors, respectively, and 
C = C L + C R . It has been shown before^ 7 ' 20 ) that C = 1 is a critical line 
for the mirror model. Only for Cr = Cl = 1/2, can the mirror model and 
the rotator model be mapped onto each other. That is, by properly replacing 
right mirrors and left mirrors with either left rotators or right rotators, we 
can produce the same trajectory with the same probability. For Cr ^ Cl, 
the trajectories can still be mapped onto each other, but the probability to 
generate the same trajectory by mirrors is no more the same as by rotators, so 
that these two models cannot be mapped onto each other anymore. However, 
we will show later in this section, that the mirror model can still be related to 
a rotator model, even when Cr ^ Cl- 

The mirror model can, for all C = 1, still be mapped onto an isotropic bond 
percolation problem at criticality, (except in the two limits: Cr = 1, Cl = 0, 
and Cr = 0, Cl = 1, where the moving particle simply zig-zags to infinity). 
Therefore the trajectory length distribution, characterized by r = 15/7 and 
the fractal dimension df = 7/4 are still the same as before. However a cannot 
be defined anymore for the mirror model at C — 1, because we are on a critical 
line rather than at a critical point, as was the case for the rotator model. 

Unlike in the rotator model, in the mirror model right mirrors can either 
be on the outer or on the inner sides of a trajectory, therefore the method 
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of separating trajectories into clockwise and counter-clockwise ones will not 
break the symmetry between the number of right mirrors M R and the number 
of left mirrors Ml on the trajectories. Indeed, our numerical results show 
that the distribution of M R /N (JV = M R + M L ) is a Gaussian centered at C R 
rather than a double Gaussian as we found in the rotator model, Fig. 13. The 
following power law describes how this distribution narrows as N increases 
(Figs 13,14), 

/ \ Mr r I ^ / \ Ml C I \ jvr-0.50 

< l^y- - > = < - Cil >~ N (16) 

The exponent 0.50 is not surprising here. The key point is that C = 1 is a 
critical line, so that for a large trajectory, we cannot tell at which concentration 
it is generated. In other words, the number M R /N can fluctuate freely about 
its mean value C R . According to the central limit theorem, we get then the 
observed result. 

Memory effects can be studied in the same way as for the rotator model, 
because each site on the trajectories can be only visited either once or twice, 
the same as for the fully occupied rotator model. Our numerical results show 
that the scaling behavior of Ni/N and N 2 /N is the same as that for the rotator 
model, except that the asymptotic value now depends on the concentration 
C R , (Fig. 15), 

< § - a(C R ) >=< b(C R ) - § >~ iV" - 57 (17) 

From Ni + N 2 = N, we have a(C R )+b(C R ) = 1. For C R > 1/2, our simulations 
show that, a(C R ) is monotonically decreasing with C R , while for C R < 1/2, 
a{C R ) is monotonically increasing with C R . At C R = Cl = 1/2, we recover 
the result of the rotator model, i.e. a(l/2) = b(l/2) = 1/2. In fact, a(C R ) can 
be fitted very well to a parabola (Fig. 16 ) except at C R = and C R — 1, 
where Ni/N — 1 and N 2 /N = indicating a discontinuous transition at these 
two limits. 
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The scaling behavior for S/N can be derived from the relation N1+N2 = N, 
Nx + 2N 2 = S, and Eq. (^), 

< I - [1 + b(Cn)) >~ N~ - 57 (18) 

This scaling behavior is the same as that found for the fully occupied rotator 
model, Eq. (|9|), except that the asymptotic value, is now a function of Cr, 
rather than a single value 3/2, as for the rotator model. 

As we mentioned before, each mirror on a given trajectory can be replaced 
by an appropriately chosen rotator to reproduce the same trajectory. After 
this operation, one can see that all right rotators are either on the outer side 
or on the inner side of the trajectory depending on the velocity direction of the 
particle, the same feature we found for the rotator model on a fully occupied 
lattice. We will refer to this rotator model generated from the mirror model 
as a quasi-rotator model. It is natural to study, in the same way as for the 
rotator model, the scaling behavior of the number of right rotators Nr and 
left rotators Nl on the trajectories for the quasi-rotator model. Our numerical 
results indicate the same power law behavior is obtain ed as in the rotator 
model, even when Cr 7^ Cl where the rotator model and the quasi-rotator (or 
mirror model) cannot be mapped onto each other, (Fig. 17). 

JV» _ 1 N, 1 ^ , 7 

N 2 N 2 y ' 

The asymptotic values for N R /N and N L /N are equal and the same as for the 
rotator model on the fully occupied lattice. This is not surprising, because 
Nr is associated with the outer side of a trajectory, while Nl is associated 
with the inner side of a trajectory (Fig. 3b). Although the mirror model for 
all C = 1 corresponds to a quasi-rotator model, it is not clear how the quasi- 
rotator model can be generated directly, other than deriving it from a mirror 
model. 
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3.2 The partially occupied lattice 

Here we only give a very brief discussion of the mirror model on a partially 
occupied square lattice. Again, this model can not be mapped onto a percola- 
tion model due to the fact that trajectories can cross themselves (Fig 2b). At 
Cr = Cl = Ce = 1/3, the mirror model can be mapped onto a growing self- 
avoiding walk trail introduced by Lyklema^ 21 ) some time ago. Previous studies 
have shown that the distribution of the size of closed trajectories and the frac- 
tal dimension of large trajectories are drastically changed as soon as C < 1: 
r = 1 and df = 2, with logarithmic corrections^ 4 ' 7 ' 22 ' 23 ^, indicating a first or- 
der phase transition. It has also been shown by numerical simulations that 
the moving particle exhibits super- diffusive behavior^ 7 ' 23 ), i.e. the diffusion 
constant increases with time logarithmically, suggesting that the hyperscaling 
relation r — I — 2/df needs logarithmic corrections. It is not difficult to show 
that no "reflector" can be formed in this model and consequently, that each 
site on a trajectory can only belong to either Ni or N 2 . 

4 Critical behavior of the rotator and the mir- 
ror model on the triangular lattice 

4.1 The fully occupied lattice 

In order to investigate the universality of the trajectory scaling behavior found 
on the square lattice, we also studied the triangular lattice. Similar to the 
rotator and mirror model on the square lattice, fixed rotators or fixed mirrors 
are randomly placed on the sites of the triangular lattice. The particle moves 
in six directions along the bonds of the lattice and its velocity is rotated by 
scatterers over an angle of ±27r/3 upon each collision, (Fig. 18). Since both 
the mirror model and the rotator model are discussed in this section, in order 
to avoid confusion , we define C r ^ tator and C r L otator as the concentration for right 
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rotators and left rotators, respectively, while we define C R lvror and C r £ irror as 
the concentration for right mirrors and left mirrors, respectively. Similarly as 
before, C represents the total concentration of the scatterers. 

For (7=1, earlier papers have shown^ 3 ' 7 ) that the mirror model and the 
rotator model can be mapped into each other by replacing all the right (left) 
mirrors on the trajectories by either right (left) rotators or left (right) rotators, 
respectively, provided C r R otator = C n R lirror and C r L otator = C™ rror . Note that this 
is a global transformation rather than a local transformation, as was the case 
with the mapping of the rotator model onto the mirror model on the square 
lattice for C R = C L = 1/2 (Fig. 19). In other words, the mirror model and 
the rotator model are equivalent . Again it can be shown that the trajectories 
cannot cross themselves and that all right rotators are either on the outer side 
or the inner side of the trajectories, depending on the direction in which the 
particle trajectory closes. 

Both the mirror model and the rotator model can be mapped onto a site 
percolation rather than a bond percolation problem, as was the case for the 
rotator and the mirror models on the square lattice. For site percolation on the 
triangular lattice, the critical probability p c to occupied a lattice site is 1/2, so 
that the critical point occurs for C^ irror = C% irror = C r L otator = C r R otator = 1/2. 
Here, we will only consider the rotator model, since the mirror model can be 
done in the same fashion. The number of right rotators N R and left rotators 
N L contained in the closed trajectories, exhibit the same power law behavior 
as on the square lattice, (Fig. 20), 

< *5 _ I > c= < I _ ^£ >„„ N -°« (20) 
N 2 2 TV v 1 

where N = N R + N L . The average is taken over all trajectories which contain 

N sites. 

For the triangular lattice, the maximum number of times a site on a tra- 
jectory can be visited by the moving particle is three. This can be shown 
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by considering again the winding angle W. Upon each collision W will either 
increase or decrease by 27r/3. Therefore, when the particle returns to the same 
site, the winding angle will be a multiple of 2n/3. However, there are only 
three allowed angles of the form 2mr/3 on the triangular lattice, where n can 
be 0, 1 or 2, so that the maximum number of times a site be visited by the 
moving particle is indeed three. We define, N±, N 2 and N 3 as the number of 
sites visited by the particle once, twice, three times and four times, respec- 
tively, while S is the length of the trajectory. These quantities satisfy the 
following two sum rules, 

N 1 + N 2 + N 3 = N (21) 

and 

N\ + 2N 2 + 3N 3 = S (22) 

The scaling behavior of Nx/N, N 2 /N and N 3 /N is very similar to that found 
on the square lattice (Fig. 21), 



< ^ - Kt >~< K 2 - ^ >~< K 3 - ^ >~ AT" - 57 (23) 
where K x = 0.3197, K 2 = 0.4052, K 3 = 0.2751. The sum of K x , K 2 and K 3 is 



one as required by Eq. fl21|) . Although all values of Ki (i = 1, 2, 3) are close to 
1/3, they are not equal. Since each site has six adjacent bonds, the sites that 
belong to Ni and N 2 form the "surface" of the trajectories, while the other 
sites that belong to A^3 form the "bulk" of the trajectories. From Eq. ([23|), we 
see that the "surface" is much bigger than the "bulk" . The asymptotic scaling 
behavior for S/N follows from Eq. (0), Eq. © and Eq. (0), to be, 

<S/N- 1.9554 >~jV-°- 57 (24) 
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4.2 The partially occupied lattice 

For the partially occupied lattice, i.e. C < 1, the rotator model and the 
mirror model can still be mapped into each other by a global transformation. 
That is, for a given trajectory, all left (right) mirrors on the trajectories can 
be replaced by either left (right) rotators or right (left) rotators respectively 
while empty sites remain empty, (Fig. 22). However, neither of the two models 
can be mapped onto a percolation problem, because the trajectories can cross 
themselves via the empty sites and consequently, the inner side and the outer 
side of large trajectories cannot be defined properly anymore, as was also the 
case for the square lattice partially occupied by scatterers. Nevertheless it 
was found numerically by Cohen and Wang( 7 ) that for C < 1, there exists 
a linear critical line Cr = Cl — C/2 (Fig. 23) and the cluster perimeter 
distribution exponent r = 15/7 and the fractal dimension df = 7/4 of the 
trajectories along this critical line belong to the same universality class as that 
of site percolation. However, the existence of this simple critical line at low 
concentrations could not be established due to the prohibitively long numerical 
calculations needed. The lowest concentration studied in our simulations was 
C R = C L = C/2 = 0.225. 

Since the trajectories can cross themselves, the winding angle W for closed 
trajectories is no more restricted to 2tt or — 2ir, so that it is natural to study 
the scaling behavior of W. Our numerical results (Fig. 24) for Cr = Cl = 
C/2 = 0.45 show the following relation, 

<|ffl|>l^ (25) 

in contrast to the stretched exponential behavior for the partially occupied 
rotator on the square lattice, Eq. fllPf) . The asymptotic value for W(N)/N is 
zero, because the moving particle has equal probability to turn to its right or 
to its left. 
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For Nr, Nl and Ne, our numerical simulations yield a different scaling 
exponent from that for the fully occupied lattice , (Fig. 25), 

< l^p - C R \ >~< |§ - C L \ >~< |^ - C E \ >~ AT-o.50 (26) 

This scaling behavior is the same as that of the partially occupied square 
lattice, Eq. flTTI). The exponent —0.50 can be understood in the same way as 
for the mirror model for C = 1. 

Using the same kind of argument, as we used above for the fully occu- 
pied lattice, i.e. considering the change of the winding angle, we can show 
that the maximum number that a site on a trajectory can be visited by the 
moving particle is three. Unlike the rotator model on the partially occupied 
square lattice, Eq. (|TJ), we find, that NJN, N 2 /N and N 3 /N now follow the 
same scaling behavior along the critical line as the rotator model on the fully 
occupied lattice, (Fig. 26), 

< ^ - Ki >~< K 2 - ^ >~< K, - § >~ (27) 

where K\, K 2 , K 3 are functions of the concentration, which sum to 1, while 
the exponent —0.57 is independent of the concentration. Our numerical sim- 
ulations also show that K\ is a not a monotonic function of C while K 2 is an 
increasing function of C and K 3 is a decreasing function of C. However, all 
Ki (i — 1, 2, 3) appear to converge to three different constants as C decreases, 
(Fig. 27). 

5 Conclusion 

In this paper, we have given, a more detailed numerical analysis than before 
of the nature of the trajectories generated in a Lorentz lattice gas model on 
both the square lattice and the triangular lattice at criticality. Our study has 
shown that the "structural" properties of these trajectories are highly non- 
trivial and yield some new universal scaling behavior probably related to the 
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scaling function exponent a, in addition to that associated with the critical 
exponents r and df found before. We hope that this study can lead to a more 
complete understanding of the complicated trajectories generated by a particle 
moving through randomly distributed scatterers on a lattice. The results are 
summerized in table I, table II and table III. 

Our study only answered a few questions, but raised many more questions, 
as yet unanswered, among which we note the following. 

1. For C — 1, why are the asymptotic values of N R and N L as well as of 
Ni and N 2 on the fully occupied square lattice equal to N/2, i.e. why is 
K x — K 2 — 1/2? What determines the asymptotic values for K x , K 2 and K 3 
on the fully occupied triangular lattice? 

2. Our numerical simulations indicate that the ( % can be as large as 4 
for the rotator model on the partially occupied square lattice and 3 for the 
triangular lattice) approach constants, as C is decreased along the critical line; 
in other words, the structure of the large trajectories hardly changes as the 
concentration of the scatterers is reduced, but why? 

3. Although we have given an argument why the exponents describing the 
fluctuations of Nr/N along the critical line, C — 1 for the mirror model, are 
all equal to 0.50, the more fundamental exponent 0.57 still remains elusive. 
We conjecture that the exact value is 1 — a = 4/7, but it is not clear how to 
obtain this from theory. 

4. Why does the winding angle W on the rotator model on the partially 
occupied square lattice exhibit a stretched exponential law behavior, while 
it exhibits a simple power law behavior on the partially occupied triangular 
lattice? 

5. We find, that the N/N approach their asymptotic values Ki with a power 
law iV~ a57 in all cases, except for the rotator model on the partially occupied 
square lattice where the exponent has a significantly different value, —0.39. 
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What is the origin of this difference? In the next paper, dealing with scaling 
behavior near criticality, we will show that a has the same value 3/7 in all cases 
except again for the partially occupied rotator model on the square lattice. 
This may indicate that the scaling behavior of the Ni/N and the exponent a 
are related to each other. 

6. How can one understand the difference between the phase diagram of the 
partially occupied rotator model on the square lattice and that on the trian- 
gular lattice? Do the "reflectors" play an important role here? 

7. How are the lattice sites visited by a moving particle, once, twice, etc. 
correlated to each other on the lattice, with respect to both space and time? 
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Figure Captions 

Fig. 1. Scattering rules for the rotator model and the mirror model on the 
square lattice (a) right rotators; (b) left rotators; (c) right mirrors; (d) left 
mirrors. 

Fig. 2. (a) a typical closed trajectory generated by a moving particle on the 
square lattice partially occupied by rotators. Here W = 6ir, N R = 10, N L = 2, 
N E = 7, Ni = 18, N 2 = 1, N 3 = 0, N 4 = 0, N = 19 and S = 20. (b) The 
same trajectory generated by a moving particle on the square lattice partially 
occupied by mirrors. Here M R = 8 and M L = 4. 

Fig. 3. (a) A typical counter-clockwise closed trajectory on the square lattice. 
Note that all left rotators and all right rotators are on the outer side and 
inner side of the trajectory, respectively. Here W = —2tt, N r = 5, Nl = 10, 
Ni = 14, N 2 — 1, N — 15 and S = 16. (b) The same trajectory generated 
by mirrors on the square lattice. Mirrors (fat solid lines) on the outer side 
of the trajectory form a bond percolation cluster perimeter on one of the two 
sublattices (dashed lines), while mirrors on the inner side of the trajectory 
form a bond percolation cluster perimeter on the other sublattice (dotted 
lines). Here M R = 6 and M L = 9. 

Fig. 4. The probability density P(N R /N) vs N R /N for the rotator model. The 
data were obtained from trajectories consisting of 1300 to 1500 sites (o) and 
2000 to 2300 sites (+). (a) both clockwise and counter-clockwise trajectories 
are included in the data. The solid and the dashed curves are described by 
the double Gaussians; 23.7 e -™58(N R /N-o.5i82f + 23 _ 7e -7058(Ar K /iv + o.5i82)2 and 

28.0 e -9856(JWJV-0.5145)a + 2 g_ ^9856(^/^+0.5145)^ respectively . ( b ) only the 

clockwise trajectories are included in the data. The solid and the dashed curves 
are described by the Gaussians, 23.7 e -7058(AW7v-o.5i82) 2 an d 2 8.0 e ~^56(N R /N-o.5U5)^ 

respectively. For increasing values of iV the curves become taller and narrower. 
Fig. 5. Scaling behavior of N R and Nl for clockwise and counter-clockwise 
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closed trajectories, respectively, on the square lattice fully occupied by ro- 
tators, C R = C L = 1/2: log 2 < N R /N - 1/2 > c vs log 2 S (o), log 2 < 
Nl/N — 1/2 > cc vs log 2 S (+). The slope of the lines through both sets of 
points is —0.57. 

Fig. 6. Scaling behavior of Aq and N 2 for trajectories on the square lattice 
fully occupied by rotators, Cr = Cl = 1/2: log 2 < Ni/N — 1/2 > vs log 2 S 
(o). Also plotted is log 2 < 1/2 — N 2 /N> vs log 2 S (+), which have the same 
values. The slope of the lines through both sets of points is —0.57. 
Fig. 7. Schematic phase diagram for the rotator model on the square lattice. 
The inner solid line extends to C = 0.65 and is tangent to the line C — 1, the 
dashed line is conjectured. 

Fig. 8. Scaling behavior of the winding W for trajectories on the square 
lattice partially occupied by rotators, for a typical concentration C = 0.9 with 
C R = 0.477, C L = 0.423. < \W/N - (C R - C L )\ > vs \og 2 S (o). The solid 
curve is described by (C R - C L ) + 7.5 x 2(- 2 - 805 ° 18 ). 

Fig. 9. Scaling behavior of N R , N L and N E for trajectories on the square lattice 
partially occupied by rotators, for C = 0.9 with Cr = 0.477, Cl = 0.423, 
C E = 0.1: log 2 < \N R /N- C r \ > vs log 2 S (o), log 2 < \N L /N- C L \ > vs log 2 S 
(+), log 2 < \N E /N - C E \ > vs log 2 S (□). The slope of the lines through the 
points is —0.50. 

Fig. 10. The two smallest "reflectors", each of six sites with one empty site, 
forming a closed trajectory on the square lattice partially occupied by rotators. 
Note that neither of the two right rotators can be replaced by either a right 
mirror or a left mirror, so the "reflector" cannot be formed in the mirror 
model. We also note that this reflector is even smaller than that reported 
previously*- 9 ' 10 -*. 

Fig. 11. Scaling behavior of Aq, N 2 , N 3 and for trajectories on the square 
lattice partially occupied by rotators, for C = 0.9 with Cr = 0.477, C L = 
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0.423, C E = 0.1 : log 2 < Ni/N - K x > vs log 2 5 (o), log 2 < N 2 /N - K 2 > vs 
log 2 5 (+), log 2 < K 3 - N 3 /N > vs log 2 S (□), log 2 < K 4 - N 4 /N > vs log 2 S 
(x), where the values of K u K 2 , K 3 and K 4 are 0.3047, 0.2671, 0.2223 and 
0.2056, respectively. The slope of the lines through the points is —0.39. 
Fig. 12. Values of K 1 (o), K 2 (+), K 3 (□) and K 4 (x) for trajectories on 
the square lattice partially occupied by rotators as a function of C. The lines 
through the points are drawn to guide the eye. 

Fig. 13. The probability density P(M R /N) vs M R /N for the mirror model 
on the fully occupied square lattice. The data were obtained from trajectories 
consisting of 1300 to 1500 sites (o) and 2000 to 2300 sites (+). The solid 
and the dashed curves are described by the Gaussians; 3i e ~ 2950 ( M n/ N ~ 0A ) 2 an ,i 

39e -4950(M K /Af-0.4)2 ^ respectively . 

Fig. 14. Scaling behavior of Mr for trajectories on the square lattice fully 
occupied by mirrors at different Cr : log 2 < \M R /N— Cr\ > vs log 2 S, Cr = 0.1 
(o), C R = 0.2 (+), C R = 0.3 (□), C R = 0.4 (x). The slope of the lines through 
the points is —0.50. 

Fig. 15. Scaling behavior of Ni for trajectories on the square lattice by mirrors 
at different Cr along the critical line C — 1 : log 2 < Ni/N — <x(Cr) > vs log 2 S : 
C R = 0.1 (o), C R = 0.2 (+), C R = 0.3 (□), C R = 0.4 (x). The slope of the 
lines through the points is —0.50. 

Fig. 16. For the mirror model on the fully occupied square lattice, a(C^) vs 
Cr (o), 6(Cr) vs Cr (+). The solid line and the dashed line are described by 
0.5 + 0.59(0.25 - C R C L ) and 0.5 - 0.59(0.25 - C R C L ), respectively. 
Fig. 17. Scaling behavior of N R for trajectories on the square lattice fully 
occupied by rotators mapped from mirrors, i.e. for the quasi-rotator model, 
at different C R : log 2 <N R /N - 1/2 > c vs log 2 S, C R = 0.4 (o), C R = 0.3 (+), 
C R = 0.2 (□), Cr = 0.1 (x). The slope of the lines through the points is 
-0.57. 
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Fig. 18. Scattering rules for the mirror and the rotator models on the trian- 
gular lattice, (a) the six velocity directions; (b) rotator scatterers; (c) mirror 
scatterers. 

Fig. 19. (a) A typical clockwise closed trajectory on the triangular lattice. All 
right rotators and all left rotators are on the outer side and on the inner side 
of the trajectory respectively and form two site percolation cluster perimeters 
on the triangular lattice. Here N\ = 9, N2 = 3, N 3 = 6, Nr = 13, Nl = 5, 
N = 18, and S = 33. (b) the same trajectory generated by mirrors on the 
same lattice. Note all right rotators have been replaced by right mirrors while 
all left rotators have been replaced by left mirrors. Here M R = N R = 13 and 
M L = N L = 5. 

Fig. 20. Scaling behavior of Nr and Nl for trajectories on the triangular 
lattice fully occupied by rotators, Cr = Cl = 1/2: log 2 <N R /N — 1/2 > c vs 
log 2 S (o), log 2 < N L /N- 1/2 > cc vs log 2 S (+). The slope of the lines through 
the points is —0.57. 

Fig. 21. Scaling behavior of Ni, N 2 and N 3 for trajectories on the triangular 
lattice fully occupied by rotators, Cr = Cl = 1/2: log 2 < N ± /N — K\ > 
vs log 2 S (o), log 2 < K 2 - N 2 /N > vs log 2 S (+). log 2 < K 3 - N 3 /N > vs 
log 2 S (□), where the values of K u K 2 and K 3 are 0.3197, 0.4052 and 0.2751 
respectively. The slope of the fitting lines are —0.57. 

Fig. 22. (a) A typical closed trajectory generated by a moving particle on 
the triangular lattice partially occupied by rotators. Note that the trajectory 
crosses itself. Here W = 0, N R = 3, N L = 3, N E = 5, JVi = 10, N 2 = 1, 
A^ 3 = 0, iV = 11 a nd S = 12. (b) The same trajectory generated by a moving 
particle on the triangular lattice partially occupied by mirrors. Right rotators 
are replaced by right mirrors while left rotators are replaced by left mirrors. 
Here Mr = Nr = 3 and M L = N L = 3. 

Fig. 23. Phase diagram for both the mirror and the rotator models on the 
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triangular lattice. The solid critical line has been computed for C > 0.6, the 
dashed critical line is conjectured. 

Fig. 24. Scaling behavior of W for trajectories on the triangular lattice 
partially occupied by rotators, for a typical concentration C = 0.8 with 
C R = C L = 0.4: < \W/N\ > vs log 21 S (o). The slope of the line through 
the points is —0.50. 

Fig. 25. Scaling behavior of Nr, Nl and Ne for trajectories on the triangular 
lattice partially occupied by rotators, for C = 0.8 with C R = C L = 0.4 and 
C E = 0.2: log 2 < \N R /N - C R \ > vs log 2 S (o), log 2 < \N L /N -C L \>vs log 2 S 
(+), log 2 < \N E /N - C E \ > vs log 2 S (□). The slope of the lines through the 
points is —0.50. 

Fig. 26. Scaling behavior of Ni, N 2 and N 3 for trajectories on the triangular 
lattice partially occupied by rotators, for C = 0.8 with Cr = Cl = 0.4 and 
C E = 0.2: log 2 < Ni/N - K x > vs log 2 S (o), log 2 < K 2 - N 2 /N > vs log 2 S 
(+), log 2 < K 3 - N 3 /N > vs log 2 S (□), where the values of K u K 2 and K 3 
are 0.3356, 0.3391 and 0.3254, respectively. The slope of the lines through the 
points is —0.57. 

Fig. 27. Values of K\ (o), K 2 (+) and K 3 (□) for trajectories on the triangular 
lattice partially occupied by rotators as a function of C. The lines through 
the points are drawn to guide the eye. 
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TABLE I: Critical behavior for C = 1 



Square Lattice 


Triangular Lattice 


Rotator 


Mirror 


Rotator=mirror 


Nr l 

< ~N 2 >c ~ 

< ^£ _ I > ^ /V" - 57 

Af 2 ^ cc iV 


< 1— ~~ u r\ >_ 

< \^-c L \ >~ iv-°- 50 


Nr 1 
^ AT 2 ^ c 

< _ I > ^ )y-0.57 
^ AT 2 ^ cc iV 


ATi 1 

< — — - > = 

^ N 2 ^ 

< I _ ^ >^ iV-0-57 
^2 AT ^ JV 


< b(C R ) - § >= iv-°- 57 


< ^ -0.3197 >~ 

< 0.4052 - ^ >~ 

< 0.2751 - § >~ iV- ' 57 


< S _ 3 jy-0.57 
^ AT 2 ^ iV 


< § - 1 - 6(C fl ) >~ iV" - 57 


< | - 1.9554 >~ iV" - 57 
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TABLE II: Critical behavior for C = 0.9 



Square Lattice 


Triangular Lattice 


Rotator 


Mirror 


Rotator=mirror 


<|f -{C R -C L )\>~2-^ 




< If | >~iv- - 50 


<\^-C R \ >= 

< \^-c L \ >~ 

< \^-c E \ >~iv- - 50 


Super-diffusion 


^ \ N * r I ^- 

< \^--C E \ >~iV- - 50 


< ^ - 0.2671 >~ 

< ^ - 0.3047 >~ 

< 0.2223- ^ >~ 

< 0.2059 - § >~ N-°- 39 


(r = 1, d, = 2) 


< 5" - 0.3308 >~ 

< 0.3669 - >~ 

< 0.3023 - # >~ N~ - 57 


< § - 2.3670 >~ iV" - 39 




<§- 1.9715 >~ iV" - 57 
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TABLE III: Critical behavior for C = 0.8 



Square Lattice 


Triangular Lattice 


Rotator 


Mirror 


Rotator=mirror 


<|f -(C R -C L )\>~2-^ 018 




< If | >~iv- - 50 


<\%-c R \ >= 

< \^-c L \ >~ 

< \^-c E \ >~iv- - 50 


Super-diffusion 


^ \ N * r I ^- 

< \^--C E \ >~iV- - 50 


< ^ - 0.2680 >~ 

< ^ - 0.2760 >~ 

< 0.2365- ^ >~ 

< 0.2195 - § >~ N-°- 39 


(r = 1, d, = 2) 


< 5" - 0.3356 >~ 

< 0.3391 - >~ 

< 0.3254 - ^ >~ N~ - 57 


<§- 2.4075 >~ iV" - 39 




<§- 1.990 >~ iV-°- 57 
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